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DIFFERENCE EQUATIONS AND THEIR APPLICATIONS 


Louis A. Pipes 


I. Introduction. 

During the last fifteen years there has been a great deal of activity 
in the applications of difference equations in the solution of many prob- 
lems that arise in statistics, science, and engineering. 

The development of high-speed digital computing machinery has moti- 
vated the use of difference equations as approximations to ordinary and 
partial differential equations. [See bibliography, 1, 2, 3]. The use of dif- 
ference equations for solving partial ditferential equations was discussed 
in 1928 in a celebrated paper by Courant, Friedrichs, and Lewy. [4]. How- 
ever, the methods presented in this paper were not put into practical use 
until about 1943 as a result of the stimulus of war-time technology and 
with the assistance of the first digital automatic computers. Problems in- 
volving time-dependent fluid flows, neutron diffusion and transport, ra- 
diation flow, thermo-neuclar reactions, and problems involving the solu- 
tion of several simultaneous partial differential equations are being solved 
by the use of difference equations throughout the country. 

Besides the use of difference equations as approximations to ordi- 
nary and partial differential equations, they afford a powerful method for 
the analysis of electrical, mechanical, thermal, and other systems in 
which there is a recurrence of identical sections. The study of the be- 
havior of electric-wave filters, multistage amplifiers, magnetic ampli- 
fiers, insulator strings, continuous beams of equal span, crankshafts of 
multicylinder engines, acoustical filters, etc., is greatly facilitated by 
the use of difference equations. The usual methods for solving such sys- 
tems are generally very lengthy when the number of elements involved is 
large. The use of difference equations greatly reduces the complexity 
and labor in problems of this type. 

Il. Nature of the Calculus of Finite Differences. 

As is well-known, the most important concept of mathematical an- 
alysis is that of a function. If to a given value of z a certain value of y 
corresponds, we say that y is a function of the independent variable x 
and write symbolically yz): Two types of functions occur in the appli- 
cations of mathematical analysis to physical problems; first there are 
functions in which the variable z may take every possible value in a given 
interval, that is, the variable z is continuous. Functions of this type are 
studied in the ordinary Differential and Integral Calculus. 
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A second type of function exists in which the variable z takes only 
the given values z,, z,, In this case the variable z is discon- 
tinuous. The methods of the ordinary calculus are not in general appli- 
cable to such functions. The Calculus of Finite Differences deals es- 
pecially with functions of the second type, but it may be applied to func- 
tions of the first type as well. The origin of the Calculus of Finite Dif- 
ferences is usually ascribed to Brooks Taylor who published his, “*Me- 
thodus Incrementorum” in 1717. He pointed out that the relations between 
the terms of arithmetic and geometric progressions, which were known for 
a long time represented the simplest examples of difference equations. 

There are many analogies between the Differential and Integral Cal- 
culus and the Calculus of Finite Differences. For example, in the Cal- 
culus of Finite Differences, the first difference, denoted by A(A)F (2) is 
defined by, 


(2.1) ACA) F(z) = F(x+h) - F(x) 


In the Differential Calculus, the first derivative of a function F(x) 
is defined by the equation, 


A 


(2.2) DF(z) Limit _ Limit 
h 


> 


It can be shown that tle ath. derivative of /(z) may be expressed in 
the form, 


(2.3) D"F(z) = Limi 


The operator A(A) is called the forward difference operator. When 
there cannot be any misunderstanding about the increment A, the differ- 
ence operator A(h) is merely written as A, 


Basic Operators used in the Calculus of Finite Differences. 


The following fundamental operators are used widely in the Calculus 
of Finite Differences. 
1. The Displacement Operator E. 


(2.4) EF (a) = F(x+h) 

2. The Forward Difference Operator A. 

(2.5) AF(z) = F(a+h) -F (2) 
3. The Backward Difference Operator V . 

(2.6) VF(z) = F(x) -F(2-A) 
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4. The Central Difference Operator 5. 


5F (x) = F(2+4)-F (a4) 


5. The Averaging Operator up. 
(2.8) uF (2) = 
6. The Symbolic Representation of Taylor’s Expansion. 


(2.9) = Flash), 
dz 


Equivalent Operators. 

By using the above fundamental definitions of the basic operators 
used in the calculus of finite differences, the following equivalent opera- 
tors may be easily deduced, 


(2.10) 
(2.11) 


(2.12) 
(2.13) 
(2.14) 


It can be shown [5, 6, 7, 8] that except for a few restrictions, these 
operators can be manipulated according to the laws of ordinary algebra, 
and many relations between the basic operators may be obtained. 

For example, we may write, 


(2.15) | D = tats) by the use of (2.10) 


(2.16) V 
(2.17) = 2sinh(AD /2) 
(2.18) = cosh(hD /2) 


An interesting and useful exposition of the symbolic calculus of op- 
erators will be found in a paper by Bickley [9]. By the use of these oper- 
ators, and their relations between them, several finite difference formulae 
for interpolation, differentiation, integration, scanning, etc. may be ob- 
tained. These formulae are of great importance in numerical analysis. 

For example, if the logarithm in (2.15) is expanded in powers of A, 
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h 2 3. 


This expression may be used to compute the derivative of a function 
in terms of its forward differences. 
The Scanning Operator. S(h) 
In the theory of scanning, if we are given a function F(z), it is re- 
quired to obtain the function ¢(z) defined by, 
(2+h 72) 
(2.20) = F(2)dz 


(2—h; 2) 


If (2.20) is differentiated with respect to z, the result is, 


(2.21) = 


= 5F(z) 


Hence, 


(2.22) = F(z) = S(A)F(x) 


Where S(A) is the scanning operator defined by, 


2.94) — 2sinh(hD/2) _ hsinh(AD/2) 
D D (AD /2) 


If d(z) is known and it is desired to determine F(z), we may write, 
2.24) F(x) = S~ (A) d(z) 


In (2.24), S~ (A) is the inverse scanning operator. 


Ill. Linear Difference Equations with Constant Coefficients. 

The most important difference equations that arise in applications 
to dynamics, electrical circuits, structural problems, etc., are linear dif- 
ference equations with constant coefficients. 

These equations are of the type, 


(3.1) (a,E" jE" +a,)yl2) = V(2) 


In the usual applications of difference equations, h = 1, so that the 
operator F is defined by, 


(3.2) E* yx) = ylx+k) 
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Difference equations are classified in the same manner as differen- 
tial equations. The equation (3.1) is called a linear difference equation 
of order n with constant coefficients, provided that the coefficients a, 
are constants. If V(x) = 0, then the equation is said to be homogeneous; 
if not, it is said to be non-homogeneous. 

Let, 


(3.3) P(E) = (a,E"+a,_ 


The general solution of (3.1) may be written in the form, 


(3.4) y= y (2) + y,(2) 


If this is substituted into (3.1), the result is, 
(3.5) yp) =0 


Let, 
(3.6) P(E)y (2) = 0, 


and 


(3.7) P(E)y (x) = V(z) 


By analogy to the solution of linear differential equations with con- 
stant coefficients, y, is called the complementary function, and Vp the 
particular integral of the differential equation (3.1). 

Determination of the Complementary Function. 


To determine the complementary function y,(z) we assume a tenta- 
tive solution of the form, 


(3.8) = Ag” 


where A is an arbitrary constant and g a number to be determined so that 
(3.8) is a solution of (3.6). It is noticed that if the operator E* is applied 
to (3.8) the result is, 


(3.9) E*y (2) = EF Ag? = g*Aq® = q*y (2), 


It is therefore evident that if (3.8) is substituted into (3.6), the re- 
sult may be written in the form, 


(3.10) P(E)y (x) = P(g)y, (x) = 0 
If the trivial solution y,(z) = 0 is excluded, then we must have, 
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(3.11) P(g) = 0 


This is the characteristic equation of the differerce equation (3.1). 
The characteristic equation is an algebraic equation of the nth. degree 
in g. The roots of this equation determine the possible values of g that 
allow (3.8) to be a solution of (3.6). 

If the roots of (3.11) are the distinct numbers q,,---¢,, then the gen- 
eral solution of (3.6), may be written in the form, 


(3.12) y 2) Agi + 


where A ,, A,,-+- A, are arbitrary constants. 

This solution is valid provided the roots g; are distinct, regardless 
of whether they are real or complex. 

If the characteristic equation (3.11) has a multiple root q, of the rth. 
order, it can be shown [7] that the part of the complementary function y_, 


that corresponds to this repeated root is, ' 


(3.13) Yo lt) = (C+ Cyr? 


where the C,,’s are arbitrary constants. 


The Particular Integral, y (a). 

The differential equations most frequently encountered in practice 
are linear and homogeneous with constant coefficients. The various meth- 
ods for obtaining the particular integrals of linear equations with constant 
coefficients have their counterpart in the theory of difference equations. 
[5, 6, 7, 8). 


The method to be used in determining the particular integral Y» de- 
pends on the nature of the function V(r). To illustrate the general proce- 
dure, let the function V(x) be of the form, 


(3.14) Vir) = Ke™* 
The particular integral Vp in this case satisfies, 
(3.15) P(E)y, = Ke™* 


In this case, to obtain the particular integral assume a solution of 
the form, 


(3.16) Yy = Be™* 
where B is to be determined. It is noticed that, 


(3.17) Ae™*e™? 
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It is therefore evident that, 
(3.18) P(E)Be™* = Ae™*P(e™) 


Hence if (3.16) is substituted into (3.15) the result is, 


(3.19) Be™*P(e™) = 
Therefore if P(e”) 4 0, B is given by, 
(3.20) B= K/P(e™), 


and the particular integral of (3.15) is 


(3.21) yp = Ke™*/P(e™) 


In general, the particular integral of (3.1) may be expressed symboli- 
cally in the from, 


l 1 
3.22 ——V(z) = P~ (E)V( 
(3.22) PE) zr) V(x) 


If V(x) is a rational and integral function of z it is possible to ex-. 


pand P~‘(E) in a series of ascending powers of A by the-substitution 
FE = (1+) and operate with this series on V(z). If V(z) is a rational and 
integral function of x, the result of this operation will give a finite num- 
ber of terms. 

In some cases it is convenient to expand the operator P~'(E) into 
partial fractions. If this is done, the required particular integral may be 
obtained in terms of indefinite sums. [7]. 


IV. Applications of Difference Equations to the Analysis of a Uniform 
Electric Network. 

An interesting and typical technical problem that illustrates the use 
of difference equations, is to determine the current and potential distri- 
bution of the electric circuit shown in Figure 1, in the periodic steady 
alternating state. 


Ty, 


Figure 1. 
The circuit of Figure 1 is a cascade connection of n boxes. Each 
box contains the Pi Network of impedances Z and admittances Y/2, shown 
in Figure 2. 


237 
75 


MATHEMATICS MAGAZINE (May-June 


Circuit Inside the kth. Box 
Figure 2. 
Each box of the network of Figure 1 contains the Pi Network of Figure 2. 
This consists of a symmetrical arrangement of two shunt admittances 
Y/2 and a series impedance Z. 

If the sending end potential V, and the sending end current /, are 
specified, it is required to determine the potential and current distribution 
of the uniform network. If Kirchhoff’s first and second laws are applied 
to the network of Figure 2, the following simultaneous difference equa- 
tions are obtained : 


1 


In order to solve the simultaneous system of difference equations 
(4.1) and (4.2), it is necessary to separate the variables. If (4.2) is solved 
for Thea the result is, 


If this equation is used to eliminate /,_ , and /, in (4.1), the result is, 


(4.4) Ving + 0 


This is a homogeneous difference equation of the second order of the 
type discussed in Section III. Its solution is of the form, 


(4.5). | V, = Ag* 


where A is an arbitrary constant and g satisfies the characteristic equa- 
tion, 


(4.6) q?-(2+ZY)q+1=0 


The two roots of this equation, satisfy the equations, 
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(4.7) 
(4.8) 
Let 
(4.9) 
then 
(4.10) = (2+ZY) = 2cosh(a) 
The Propagation Constant of the Network. 
(4.11) a = cosh™ (1+ZY/2) = Propagation Constant. 


(4.12) 


0 


k+1 


The Characteristic Impedance of the Network. = Z, 


(4.11) Z, = 2tanh(a/2) 
y 


Determination of the Arbitrary Constants. 


The two arbitrary constants 4, and A, may be obtained from the 
boundary conditions. 


(4.15) 
( 1.16) 


Hence, 


(4.17) 


(4.18) A, = 


If these values of A, and A, are substituted into (4.12) and (4.13) 
the result may be expressed in the form, 


(4.19) V ,cosh(ka) Z sinh(ka) 
k=0,1,2,+n 


(4.20) I, = Sinhtka), 1,cosh(ka) 


4 


0 


V,=A,e%*+ A,e~™ 
V,=A4,+A, 
A,) 
liy 
If 
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cosh(a) Z oSinh(a) 
(4.21) {T] = = Transmission Matrix. 


sinh(a) 


Vy 
(4.21) (y), = 
I, 


Then, 
(4.22) (y), = (TI-*(y), 


cosh(a) 


and 


(4.23) (y), 


V. Longitudinal Oscillations of a Long Uniform Train. 

Another typical example of the use of difference equations is the 
application to the dynamical problem of determining the resonant fre- 
quencies and modes of oscillations ofa chain of n cars coupled by springs 
that satisfy hooke’s law as shown in the Figure 3. 
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Idealized Uniform Train of n Cars 
Figure 3. 

Figure 3 represents an idealized train of n cars of equal mass i 
coupled to each other by identical springs of spring constant &. The ef- 
fects of friction are neglected. The simplest method of analyzing the be- 
havior of this system is by the electrical analog method which makes use 
of the results of Section IV. [See 10 in the Bibliography] 

In this section, the dynamical system of Figure 3 will be analyzed 
by the use of a chain of differential equations and difference equations. 

If Newton’s second law of motion is applied to the rth. car in the 
chain of cars of Figure 3, the following differential equation is obtained, 


(5.1) ka r= 2,3,++-(n-1) 


In order to study the oscillations of the system, assume a solution 
of (5.1) of the form, 


(5.2) z, = y,Sin(wt+ 6) 
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where y, is an amplitude function and » is an angular frequency to be 
‘ determined. @ is an arbitrary constant. 

If (5.2) is substituted into (5.1) and the result divided by sin(wt+@), 
one obtains. 


(5.3) = 0 
where 
(5.4) 6 = (1-w?m/2k) 

To solve the difference equation (5.3) assume, 
(5.5) y, = Aq’ 


If this solution is substituted into (5.3) the following characteristic equa- 
tion is obtained : 


(5.6) q?-209¢+1-=0 
Let the two roots of this equation be ¢, and q, and let, 
(5.7) q,=e* 


Therefore, 
(5.8) (q ,+9,)/2 = (e%+e~%)/2 = cosh(a) = 0 


The solution of the difference equation (5.3) may be written in the 
convenient form, 


(5.9) y, = B ,cosh(ar) +B,sinh(ar), r= 2,3, 4, ++ (n-1) 


where B , and B, are arbitrary constants. 
The Boundary Conditions. 

In order to obtain the constants B, and B,, use must be made of the 
boundary conditions at the front and rear of the train. It is realized that 
there are no cars corresponding to r = 0 and to r = (n+1). In order to ex- 
tend the solution (5.9) to take into account r = 0 and r = (n+1) assume 
that there are phantom cars that correspond to these indices. Since the 
front and rear of the train is able to move freelii, there cannot be any 
strain in the phantom springs connecting the car r = 0 to r=1 or the spring 
that is supposed to connect the car r =n to the car r = (n+1). This condi- 
tion may be expressed mathematically by requiring that, 


(5.10) Yo= 9 
(5.11) 


and (5.9) can be used as if it were valid for r = 0 and r = (n+1). 


The boundary conditions (5.10) and (5.11) lead to the two equations. 
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(5.12) B = B ,cosh(a) + sinh(a) 
(5.13) B ,cosh(an) +B,sinh(an) = B ,cosh[{a(n+1)} +B ,sinhla(n+1)] 
B, can be expressed in terms of B , by (5.12) in the form, 


(5.14) B, = -tanh(a/2)B , 


If this expression for B, is substituted into (5.13) and some algebraic 
reductions are made, the result is, 


(5.15) B ,tanh(a/2)sinh(an) =0 


For a non-trivial solution, we must have B , = 0. Hence the possible 
values of a must satisfy the transcendental equation, 


(5.16) tanh(a/2)sinh(an) = 0 


This equation is satisfied provided that, 
(5.17) sinh(an) = 0, 


or 
(5.18) jms/n, 8 j= (-1)* 


If s is given the value s = n, (5.16) leads to an indeterminate form 
and it is not satisfied. 
The Resonant or Natural Frequencies of the Train. 

To determine the possible natura] angular frequencies of the system, 
we may use (5.8) written in the form, 


(5.19) cosh(a) = cosh(a,) = 6 = (1-w?m/2k) = cosh(jzs/n) 


Hence, 
(5.20) w3m/2k = 1-cosh(jms/n) = 1-cos(7s/n) = 


Therefore, 
(5.21) w, = 2(k/m)* sin(ws/2n), s = 0, 1, 2, 3, (n—1) 


These are the n natural angular frequencies of the principal modes 
of oscillation of the train. The frequency w, = 0, corresponds to a rigid 
body translation of the entire train. 

The Natural Modes of Oscillation of the Train. 

To determine the natural modes of oscillation of the system, write 

(5.9) in the form, 


(5.22) y, = B ,cosh(ar) +B ,sinh(ar) 
= B [cosh(ar)-tanh(a/2)sinh(ar)], r= 1,2, 
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To each value of a = a, = jas/n, there corresponds an amplitude func- 
tion, 


(5.23) ¥p, = 8), 
where 
(5.24) A(r, 8) = 


r =1,2,3,+-n (cars) 
s = 0,1, (modes) 


The general solution of the chain of differential equations (5.1) with 
the appropriate boundary conditions can now be obtained by substituting 
(5.23) into (5.2) and summing over all the possible modes s. 

The general solution of the chain of equations (5.1), is therefore, 


s=(n-1) 
(5.25) B ,t+0,) 
n 2n n 


s=0 


The index r refers to the car number and the index s pertains to the 
mode. The quantities 8, and @, are the 2n arbitrary constants of the gen- 
eral solution. They may be obtained if the initial displacements z,(0) and 
the initial velocities 2/(0) of each car is given. 


VI. Difference Equation Method for Solving Partial Differential and Ordi- 
nary Differential Equations. 

There exists a vast literature in which the solution of ordinary and 
differential equations may be solved by the device of transforming the 
differential equation into an approximate difference equation. (1, 2, 3, 11] 

From the definition of the derivative operator D and the operators A, 
V, and 4, it is easily seen that the following operators are approximations 
for the derivative operator D, provided A is sufficiently small. As A+0 the 
right members of these equations tend to D. 


:A 
(6.2) D= 
(6.3) Dez + 


' To illustrate the use of these approximations to the operator D, let 
it be required to solve the differential equation, 
2 
(6.4) =2, or D*y=z. 
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Let the initial conditions be 


(6.5) y(0)=0, and Dy0)=1 
For a sufficiently small h, we have approximately, 


Hence the differential equation (6.4) may be replaced by the approxi- 
mate difference equation, 


h? 
The initial conditions can be replaced by, 
(6.8) = 0, 


Let the difference equation be required to be satisfied at the successive 
points, z, = h, 2, = 2h, ---.2, = kh. Then the difference equation (6.7) 
may be written in the form, 


(6.9) Yew 1 + (hk = 1,2, +++) 
where 
(6.10) y, = = ylkA). 
The initial conditions take the form, 
(6.11) Yo =9 =A 
The solution of (6.9) is 
(6.12) y, = (A-h®/6)k + h®k*/6 
This corresponds to, 
(6.13) = a, 


The exact solution is, 


(6.14) yz) = 2+ 2°/6 


If A is so small that h? can be neglected, (6.13) tends to (6.14). 
Partial Differential Equations. 

The Poisson differential equation will be taken as an example of 
the use of difference equations to solve partial differential equations. 
[3]. This expation has the form, 


(6.15) (Di+D y) = fiz, y) 
In this case we have approximately, 
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(6.16) D2 = (E,-2+E>')/h? 


2 —1) 7,2 
(6.17) Dy = 


where, 


(6.18) E ula, y) = ula+h, y), E julz, y) = ule, y+h) 


If the operators (6.16) and (6.17) are used, the partial differential 
equation (6.15) may be written in the approximate form, 


(6.19) wlar+h,y) ulz, - 2ulz, y) + ulz,y-h) = 


The solution of this equation with prescribed boundary conditions is 
solved by Rosenbloom [3}. 


VII. Additional Applications. 

Additional applications of Difference equations have been made to 
the study of intermittent control systems [12] and the analysis of systems 
with hereditary characteristics [13], as well as in certain problems in the 
theory of probability 
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Editor, MATHEMATICS MAGAZINE, — 

Success to Professor D. H. Hyers, and praise for his establishment of 
the S.I.M.N.T.! (See Mathematics Magazine, -March-April, 1959.) Please 
enroll me as a member, and for my Habilit&tationvorschlag I propose to 
replace “denumerable” and “enumerable” (which sounds too much like 
“innumerable”) by “countable”. 

Let the Executive Committee also take under consideration words 
whose primary meanings do not suggest their mathematical significance, 
such as “field”, “ring”, “loop”, “histogram”, “regression”, etc. 

And can we get rid of the word “imaginary”, so misleading to-day, 
and substitute the single letter, “:” (or “j”), both as noun and as adjec- 
tive, with the definition i? = -1, and drop the V-1 ? 

William R. Ransom 
Tufts University 


EULERIAN NUMBERS AND POLYNOMIALS 


L. Carlitz 


1. Introduction. Following Euler [6, pp. 487-491], we may put 


_ 
(1.1) (A#1), 
n=0 


where H_, = H, {A} is a rational function of A; indeed 


(1.2) R LA) (A-1)"H_[A] 
is a polynomial in A of degree n-1 with integral coefficients. If we put 


n 


s=] 


then the first few values of A, are given by the following table, where 
n denotes the row and s the column; 


ll 

26 

57 57 


Alternatively, Worpitzky [15] showed that the A, . may be defined by 
means of 


n 
(1.5) 


The rational functions 1, were studied in great detail by Frobenius 

[7], who was particularly interested in their relationship to the Bernoulli 

numbers. More recently Vandiver [14] has also made use of this relation- 

ship to obtain new properties of the Bernoulli numbers. Other recent oc- 

currences are [1], [2], [12]; generalizations occur in [4], [5], [13]. In view 
247 


a 
(1.3) R (n> 1), 
= 
l 
(1.4) 
l 
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of the long history of H,, and A, it is rather curious that, on the whole, 
these quantities are not very well known. Indeed an examination of Math- 
ematical Reviews for the past ten years will indicate that they have been 
frequently rediscovered. Actually there is no detailed discussion of H,, 
in any book. On the other hand, Riordan, in his recent book [11], does 
develop a few basic properties and indicates the connection of A, , with 
certain combinational problems. 

The present paper is mainly expository. We include numerous proper- 
ties of H, and the related polynomial 


n 
(1.6) A) = (") Hu", 
r-0 


indicate the connection with Bernoulli numbers and polynomials and final- 
ly obtain some arithmetic properties of H,. For the combinatorial appli- 


cations the reader is referred to Riordan’s book [11]. The H,, also occur 


in certain criteria for Fermat’s last theorem; this is discussed at length 
in Bachmann’s book [3] and will not be considered in the present paper. 


2. The defining relation (1.1) is evidently equivalent to 


(2.1) (H+1)"=AH, (n>0), Hy=1, 


n 


where after expansion of the left member, H” is replaced by H_; we shall 


n°? 
use this convention frequently. If f(z) is an arbitrary polynomial in 2, 
(2.1) implies 


(2.2) f\H+1) = Af(H) + (1-a) (0). 
In particular, for f(z) = (7), we get 
(2.3) (m>0), 
which implies 
Repeated application of (2.3) gives 


in particular we have 


(2.5) 4) = 1. 


1959) EULERIAN NUMBERS AND POLYNOMIALS 249 


It also follows from (2.2) and (2.3) that 


(2.6) (Mery (m2); 


in particular by (2.5) 


(2.7) 
mn) (1-A)™ 


Again, if f(z) is an arbitrary polynomial of degree n, we recall that 


n 
fle) = > 


r=0 


where 


Afla) = = AT (2). 
Using (2.5) we get : 


n 
(2.8) fH) = 


r=0 
Since 


r 


ATA0) = 


we may write 


(2.9) AH) “Yon Als). 


s=0 


It is perhaps of interest to mention that (2.9) can also be obtained as 
follows from (1.1). For |z| sufficiently small we have 


e7-r 
oo n r 
Yours 8", 
n! 


n=-0 r=0 s=0 
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so that 


n 
(2.10) Hy = 8" 
r=0 s=0 


and (2.9) follows at once. 
Returning again to (1.1) it is easily verified that 


(2.11) R [A] = 

also differentiation yields 

(2.12) Ray, = (me DAR, 

where FR, is defined by (1.2). From (2.11) and (2.12) we immediately ob- 
tain 

(2.13) A 


and 


(2.14) A = 8A, ,_,+(n-s+2)A,,. 


Nn+1,8 n,s 


By means of (2.14) one can easily extend the table (1.4). A convenient 
check is furnished by the formula 


n 
= > Ane (n>1). 
s=1 


We also note that &,.[-1] = 0 for n>1. 

Frobenius remarks that it follows from (2.12) that the n-1 roots of 
R,{A] = 0 are real, negative and distinct; also for each root A, the recip- 
rocal A, is also a root. Moreover the roots of RF, [A] = 0 are separated 
by the roots of #, [A] = 0. In the next place, by (2.10) 


n 

R,, = avery (-1)"-4( 
r=0 j=0 
n r 


r=0 s=0 j=0 
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n n—-s 
Sevres 
s=0 j=0 r=} 


n-s 
s-0 j=0 
Hence, by (1.3) and (2.13), we get 


(2.15) Ang = 
j=0 
a formula due to Euler. 
Since 


n 


n 
f= 


r=0 0 


it follows from (2.10) and (2.15) 


n 


n 


n 
s=1 s=] 


which establishes (1.5). This proof is taken from Frobenius. 
3. We now consider the polynomial 


(3.1) H,(u) = A) = (u+H)” 


defined by (1.6). We evidently have the generating function 


1-A 


n=0 


It follows from (3.2) that 
(3.3) H,(u+1)-AH,(u) = (1-A)u” 


Moreover (3.3) uniquely determines the polynomial H,(u). If fiz) is an 
arbitrary polynomial in 2 then by (3.3) 


flu+1+H) -Aflu+H) = flu), 
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from which it follows that for given f(z) the difference equation 
(3.4) glutl) —Ag(u) = (1-A)flu) 
has the unique solution 
glu) = glu| A) = flu+H). 
It follows at once from (3.2) that 


(3.5) H/(u) = nH, _ ,(u) 


and generally 


= ni (n>), 


(n-r)! 


which implies 


n 
(3.6) H(utv) = 
r=0 


If we differentiate (3.2) with respect to A we get 


} 


which reduces to (2.12) for u = 0. We also note that (3.2) implies 


(3.8) = dd. 


We remark that by (2.10) and (3.1) we have 


n r 


r=0 s=0 


Again, if we put 


m—| 
(3.10) f Au, m) = 
j=0 


then it follows from (3.3) that 


H,(u+m)-AH, (u) 
1-A 
The polynomial /,(0, m) is usually called a Mirimonoff polynomial; a more 


(3.11) f,(u, m) = 
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general polynomial is discussed by Vandiver [14], see also Bachmann 
(1, p. 117). 
To get a multiplication theorem for H,(u) we consider 


00 m—\ 


n: m 
n-0 r=0 r=0 


1-A™ 


1-A n! 
n=0 


nimu | A), 


which yields 
m 1-A 
r=0 
An interesting special case of (3.12) is obtained by taking A = ¢, where 
Then (3.12) reduces to 


(3.13) m™ oH, (ur | é) =H 


r=0 


Nielsen [9, p. 54] has proved that the multiplication theorems for the 
Bernoulli and Euler polynomials characterize those polynomials. Suppose 
now that f(u|A) is a polynomial in u of degree n that satisfies the equa- 
tion 

m-1 
m 
r=0 


for some value of m>1. Put 


n 
flu| A) = > A 
s=0 


Then by (3.12) and (3.14) 
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n n 
s=0 s=0 


This requires that 
(3.15) [A™]=A,[A] (O<s<n). 


Now assume that (3.14) holds for two values of m>1, say m 


and 
Then it is clear that (3.15) becomes 


A,A™) =A, A,[A)=0 (<s<n). 


flu| A) = A, [AJH,(u| A). 


For (3.13) the situation is somewhat simpler. If g(u| 2) is a polyno- 
' mial in u satisfying the equation 


(3.16) m” (uk | ¢) = g(mu| 2), 
r=0 


where = 1, #1, then if we put 


n 
2 >) A(OH,(u| 
s=0 
and assume that (3.16) holds for one value of m>1, then it follows readily 
that 


gulO = (ul 
where A,(¢) is arbitrary. 


It may be of interest to mention also an addition theorem satisfied by 
H,(u| A). Since 


2 
1-A , eu L+A 1-A | 


it follows at once that 


n 
(3.17) (*)H(u| AH, _,(o|-A) = | 
r=0 


Therefore 
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We note also that from the identity 


e7-dr e7+A 


follows 


(3.18) (1+A)H | (uw | -A) = (¥. a?) 


while from 


e7-A e7=14A 


follows 


n 
(3.19) (aS | ADH, v| 1-A) = AH, (u+ u+v!1—A). 
r=0 


4. It is familiar that the Bernoulli polynomial B (uv) may be defined 
by [10, Chapter 2] 


(4.1) 


then 


where ¢ is a primitive m-th root of unity and m+r. This evidently implies 


r 


(4.2) ng 


For m|r, on the other hand, we have the multiplication theorem for Ba): 


(4.3) B,, = B, (mu) 
s=0 


re" 
e*-1 
n=0 
oo 
n=0 s=0 s=-0 
r 
e 
s=0 
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Multiplying both sides of (4.2) by ¢%, summing over r and using (4.3), we 
get 
H,_ ,(mu| 


(4.4) n"B ) = B,(mu) 
r=1 


where 0<t<m. 
We recall that the Bernoulli function B,(u) is defined by 


(O<u<l), B,(u+l) = 
Similarly we define 


where ¢ is some root of unity 4 1. With these definitions of H (u| ¢) and 
B, (2) it is easily verified that the formulas (3.8), (3.12), (3.13), (4.2), 
(4.3), (4.4) hold for the barred functions; in particular in (4.4) the restric- 
tion 0<¢<™m is no longer necessary. 


We remark that for m even and ¢" = -1, (4.2) reduces to the known 


formula 


m” 3 (-1)°B,, (+) = , 
s=0 


where E(u) is the Euler polynomial of degree n-1. 


For €" = w, where w7+ +1 = 0, m = 3, u = 0, (4.2) becomes 


= 3" |B, +07 


For n even, it is known that [10, p. 22] 
B,(4) = B,(2) = 
from which it follows that 


(4.6) H,,_ = (3"-1)B, (n even) . 


On the other hand for n odd > 1 we get 
(4.7) lo] = 3"wB,(4) 


Again for ¢" = i, m = 4, we get 
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It follows that 


(4.8) ni-1)H,_ = 2%2"-1)B, (n even), 


(4.9) n(i+i)H,_ [i] =-29°B,(4) (n odd > 1). 
For A = —1 we. have [10, p. 28] 


B 
(4.10) H 1) = 21-27. 


5. We now obtain some congruences satisfied by H,. If in (2.8) we 
take f(z) = 2“1-2”)", we get 


n 
s=0 
Assume that p is a prime such that 
(5.1) (p-1)p°—* | w; 
then by Fermat’s theorem 
s“(1-s”)* =0 (mod (p”, 


We thus obtain (Frobenius) 


k 
(5.2) (~1)*¥-( (mod (p", )) 
s=0 


valid for k>0, n>0, provided w satisfies (5.1). This result is referred to 
as Kummer’s congruence for H,. In (5.2) A may be an indeterminate or an 
algebraic number such that (p, 1-A) = (1). In particular A may be an /-th 
root of unity, where / # pf ‘ 


By means of (5.2) it is proved in [5] that the coefficient A, . of (1.3) 
satisfies the congruence 


(5.3) A A,, , (mod 


n+b,3 


where n>e, b = 
Another interesting result, also due to Frobenius, is 
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(5.4) H,, = =} (mod 


where again w satisfies (5.1). Indeed we can prove a slightly more gen- 
eral result by using (3.11) with n = w, u = 7, m = 0 (mod p°). We get 


m—\ 


i-0 


A-1 


where the integer r, is defined by 
(5.5) r=-f, (mod p) (0<r,<p). 
We have finally 
(5.6) = 1-A? “te (mod p*), 


which reduces to (5.4) for r = 0. 
The interesting congruence 


(5.7) = (mod p) 


is due to Mirimanoff [8]. It can be proved rapidly as follows. From (3.11) 
we get 


2 
j=0 


so that 


(5.8) = AX 


Replacing A by 1-A, we get also 


[1-a]) = 


d 


Consequently 
= 
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basic principles. 
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Reference book for anyone working in mathematics or its 
applications. 


A text book for survey (general) college courses in mathematics. 


THE FIRST SEVEN CHAPTERS 
A cure for those who have the mistaken notion that they “just 
cannot get mathematics”. 


Preparatory reading for students planning to specialize in mathe- 
matics, physics, engineering or other 


A text book for general courses in school (might be time 
for additional exercises). 7” * 


Library of Congress number 58-8471 


From the Preface to: 
“THE TREE OF MATHEMATICS” 


There is a great deal being published about mathematics these 
days, and that is fine; but this book is mathematics in the sense 
that it presents the epitomes of the main branches of the subject 
beginning with high school algebra and extending far into graduate 
work. 

When writing this treatise the authors gave great attention to 
making it both meaningful and understandable, an art in which most 
of them are pastmasters. The practice was to start on ground familiar 
to everyone and construct a highway, free from road blocks, through 
the wonderful world of mathematics. 
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where C is independent of A. Clearly C = 0 and (5.7) follows at once. 
Again by (3.11) 


2 
[Al = S (j+ 

j=0 
so that 


d 2 
Then by (5.7) 


(AP-A)R = AP(1-A)R 


; 


now using (5.8) and (5.9), we get 
-5A%1-a) 2R = +A7(1-a)?*'R 


so that 


where of course p> 2. This result also is due to Mirimanoff. 
We note that from (2.1) and 


(?=") = (-1)" (mod p) 


it follows that 


p-l 
(5.11) = AH,_, (mod p), 
r=0 


p-l 
1 
r=1 


Finally from (3.11) follows 


(5.13) = (1-)FR, (mod p), 


r= 


where p = 2k+1 and (r/p) is the Legendre symbol. 
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PROPORTIONAL METRICS IN V VARIABLES 
H. Randolph Pyle 


Introduction. 

Suppose that we are given a surface S with metric ds? = Edx*+ 
2Fdzrdy+Gdj? and a surface S’ with metric ds’? = E ‘de’? +2F dx ‘dy’ + 
G‘dy’*. The point (x’,y’) of S’-is related to the point (z, y) of S by the 
transformation 2’ = y), y’ = y). The necessary and sufficient 
condition for the conformal mapping of S on S’ is ds’? = A*ds*. We have 


shown in a previous paper (1) that this is equivalent to a set of linear 
relations between the partial derivatives of ¢ and w that reduce to the 
Cauchy-Riemann equations for isometric parameters. It is the purpose of 
this paper to find similar relations when the metrics are expressed in 
terms of n parameters instead of two. This will be accomplished by the 
use of matrix methods. 


Conditions for proportional metrics. 


We shall express the first metric as ds? = a; and the 


second as ds”? = js and impose the condition that ds” = 


\*ds*. The y’s are related to the z’s by the transformation ¥;, = Oley, 
z,), (¢ = 1,-+n). We shall make the usual assumptions concerning the 
existence and the continuity of the partial derivatives of the ¢;, and that 


the Jacobian does not vanish. Let the symbol ¢;, stand for at, so that 
z 
k 
. The Jacobian J is the determinant of the matrix 


$11 Pig Pin 
A=[ %2, Pan 


Pan 


Or, Ox, 
(Al 
Png 
Let 
261 
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(a;; = Aix bis = b) be the matrices of the coefficients of the differentials 
inds* and ds “* respectively. Designate the determinants of these matrices 
by {|A|, |B|, and assume that |A| # 0, |B| 4 0. Let dX = (dz ,, dz,,++-dz,), 
dY = (dy,,dy,,++ dy,) be row vectors and (dX), (dY)™ be the corres- 
ponding column vectors. 


Using matrix multiplication, we have ds? = (dX)A(dX)" and ds® = 


(dY)B(dY)™. Since dy; = = MGX)™ and = (dX)A™. 
=1 


Then ds“? = (dX) \7BMdX)*. When ds” = ATBA = A? A, 
Let $i; represent the cofactor of $3; in J, and 


n 
is the identity matrix. Now let L;; = L 
an 


Since ATBA = 472A, A(ATBA) = A(A?A), or (AAT)(BA) = A2(AA), and 
JL = Also = But |AT| =|A]=J, and |A?4| =A2"|Al, 


2 
so that J?|B| = A?"\A], and A? = 


Since JL = \2L, if we call » = 
and Li; = or 


|A} 7(N=—2)/Nn 


a=? we have L = pL, 


One Pin Pon Pan 
Then 
0 wee G 
AAT =[0 J O 
was 
( 21 22 -) L=AA 
n n 
k=1 k=1 
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as the necessary conditions for ds’? = A7ds?. 
Since each of these steps is reversible when J, |A|, |B| 4 0, these 
conditions are also sufficient. 


Special cases. 
1. When n = 2. 


ds* = Eda? +2F dz 


ds”? = E ‘dy? + 2F ‘dy dy, +G ‘dy? 


EF 2 2 
>, oy Vy 


Li; = BL; ; gives us the relations previously derived by other 
methods (1). 
2. In an isometric coordinate system, aj; = bi; = 8 (Kronecker delta), so 
n—2 
that |A| = |B) =l,p-J", 


When n = 2, these become the Cauchy-Riemann equations. 


Let W = (y,,---y,) be a variable vector for which y; = ¢,(z ,, ++ z,). 


Then - W; = (¢ The outer product of the n-1 vectors 
W 


jer Ws, is a vector Q; whose components are (3 
When 


$i; 
But Q; is orthogonal to the vector space spanned by the n-1 vectors men- 
tioned above, so that W 5 is orthogonal to this space. Hence the set of 
vectors W .,W,,---W, is an orthogonal set. 
1, H. Randolph Peed Conformal Mapping of Surfaces. Duke Mathematical Journal. Vol. Il, 
Whittier College, Whittier, California 


CARRY-OVER 


Colonel John W. Ault 


It is probably safe to say that every mathematics teacher believes 
that any student will derive great benefit from a broad background in 
mathematics even though his life work may be in a field which requires 
only a basic skill in arithmetic. This feeling is based on the fact that 
there is a certain amount of “carry-over” of the mental discipline to which 
he has been exposed. 

Without any attempt to open the whole subject of carry-over, since it 
is one that most people consider debatable, the following remarks are de- 
voted to the simple approach to problem solving. By the simple approach, 
we mean the basic procedure of isolating the question to be answered and 
listing the given facts which can be brought to bear upon the question. It 
is far too evident that most students believe that (a) Given, (b) To Find, 
and (c) Solution approach to a problem should be left behind at the com- 
pletion of high school plane geometry. 

It is believed that if the student can be made to idealize this basic 
approach to problems, whether they be in mathematics, physics, econom- 
ics or whatever, that there will be some carry-over. In a military school 
it is normal to require that a student’s room be kept in apple-pie order at 
all times, but if he does not build up in his mind the fact that neatness is 
a better way of life and idealize this trait to the point where he is slight- 
ly dismayed by others who do not share this belief, then his first bachelor 
quarters after graduation are pretty apt to be sadly disorganized. 

Just how the teacher can build up this idea to the ideal level is 
hard to say. Of course, great benefit can be derived from constant em- 
phasis of the principle in lectures and drill sessions, and grasping every 
opportunity to show the student who has difficulty with a problem that he 
would have avoided his difficulty by using this approach. Perhaps in the 
long run the most effective teaching device will be the fact that the 
teacher himself applies this method to all problems met in mathematics 
or elsewhere. 


United States Air Force Academy 
Colorado 


A NOTE ON SIMPLE CORRELATION 
Maurice Fréchet 


Translation from the French 
by 
C. de la Menardiere 
Department of Modern Languages 
University of Alabama 


In Volume 22, No. 2, 1958, pp. 57-69 of Mathematics Magazine, ap- 
peared an article by C.D. Smith entitled, “On the Mathematics of Simple 
Correlation.” The aim, as given in the article, is to state the problem as 
it usually occurs and show how the measures of correlation may be de- 
rived by elementary algebraic calculations. This aim was completely ful- 
filled. The purpose of this Note is to add some comment regarding inter- 
pretation of measures of correlation which follow from theoretical con- 
siderations. 

The paper states, p. 58, “The method of correlation is designed to 
measure the strength of common influences between two variables.” We 
may note first that the coefficient of correlation r does not always indi- 
cate what influences may exist. Yule, and others, (1), have found cases 
where the measure has a value near its maximum where no common in- 
fluence is known to exist between the two variables under consideration. 

In fact, the thing that can best measure the coefficient of correlation 
is the thing that can measure a good index of correlation. We seek the 
total accuracy of the determination of one of the two variables when the 
other is given. In other words, an index of correlation has nothing to say 
as to the causes of this accuracy or the lack of accuracy, not even the 
existence of the causes. It gives a numerical fact. For instance X, is 
more or less determined by knowing X,. This is why the term, “nonsense 
correlation,” seems deceiving to us. In all cases, the value of an index 
of correlation has a meaning if the numerical data from which it has been 
extracted have been correctly taken. But that meaning is only one of the 
characteristics of the numerical table of correlation. 

When a good index is equal to 1, we may state four conditions, name- 
ly: 


I. To all values of X, there corresponds only one value of X ,. 
Il. Conversely when it is equal to 0. 
Ill. X , is constant when X, varies. 
IV. And conversely. 
Moreover it is known that the classic coefficient, r, verifies only two of 
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these conditions (I and IV) (1). In some cases it is a bad index, the use 
of which must be controlled. The value of a good index of correlation 
does not dictate to the statistician the above conclusion. It invites him 
to look for the causes of the mathematical dependence, when for instance, 
the index is close to one. 

Consider the coefficient of linear correlation given by the mathe- 
matical formula: 


No 


(1) f= 


Also consider the correlation ratio, 7 as defined by n = om, 


(2) We ask: 1. Why is it a correlation measure? 
2. Why 7 is introduced if r is already known? 
Karl Pearson, having noticed the faults above concerning r, tried to 
make a new formula. 
A calculation, quite elementary, permits one to see better the pass- 
ing from r to 7 (1). 


Indeed if one has p = ‘> one has according to (1) and (2) 


22,2, 


Then since =z, = 0, 
,-m)z, 
where m is the average of X, for z, = x, and 0, is the variance of m. 
From this it follows that: 


r= pn 

Here p is simply what becomes r when one replaces the set of points 
X,, X, by the line of average X ,, when X, is given. Then p is, like |r}, 
at the most equal to | in absolute value and one has |r| < 7 < 1. One sees 
that p does not vary when the set of points (X ,, X ,) stretches or expands 
around the line of average. . 

Consequently the coefficient r contains a strange factor, p, which 
tells us nothing about the proposed sum. Is the set spread or contracted? 
Are the values of X, more or less determined when one knows \.? By 
eliminating that strange factor, one has this time an index n, which veri- 
fies three of the conditions (I, II, II], IV). And even if it does not verify 
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exactly the condition III, which is to say, if when it equals zero, X, is 


not necessarily constant, at least the average of X , is constant, for given 
X,, when X, varies. 

Sometimes one considers that 7 is much more sensitive than r to the 
grouping of data by intervals and some uncertainty arises from grouping. 
In my opinion, this fact indicates superiority of n over r. Because the 
grouping by intervals modifies not the real correlation between X , and X, 
but that which is represented on the correlation table. An intex of corre- 
lation should reflect this change and consequently warn the statistician 
not to much reduce nor to much spread the grouping intervals. 

It. would be useful to add that, in spite of the superiority of » over 
r, neither r nor 7 will satisfy the four conditions above. Other authors, 
pursuing the improvement realized by Karl Pearson, have defined several 
correlation indexes satisfying the conditions I, II, III, IV. One such 
_measure for instance is the index called connection by Gini (2). It would 
be desirable, for practical purposes, that the conclusions drawn from the 

alculation of r and » by means of the calculation of the index of connec- 
tion of Gini be controlled. The numerical calculation was considerably 
simplified by Salvenini. 

We have explained in what manner r is inferior to n (itself insuffi- 
cient) establishing by an elementary and rapid calculation the formula 


r= 


Here p is the linear coefficient of correlation that one could obtain re- 
placing all the values of Y, when X is given, by their average. It is a 
geometric characteristic of the curve of the averages, equal to 1 when 
that curve is a straight line, and, it is only in the calculation of y that 
the dispersion of Y intervenes when X is given. 

Since the coefficient r gives satisfactory results when the couple 
(X ,,X,) satisfies the Laplace-Bravais Law, called normal to two varia- 
bles, one expects to go back to that case, by transformation of the scales 
of X, and X, separately, in order that each of the new variables satis- 
fies the, “second Law of Laplace,” called normal Law. We have shown 
(3) that even after this transformation, the intended result has not been 
reached; the new law of the couple (X ,, X,) modified is not necessarily 


the binormal law. When the frequency laws, called marginal, of X , and 


X,, considered separately, are given, the linear coefficient of correlation 
r of any tables of correlation having those margins remains included be- 
tween two numbers r,, r,, which are generally not equal either to —1 nor 
to 1, but simply included between -1 and +1 (4). These values of r,, Pes 
of contrary sign, correspond to two monotonic functional relations R me R, 
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between X, and X,. And reciprocally, if a monotonic functional relation 


exists between X, and X,, it is identical to R, or R, (5). 
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SOME FURTHER NOTES ON THE THEORY 
OF CORRELATION 


C. D. Smith 


The paper, referred to in the Note of Fréchet, proceeds from the Re- 
gression Model as given by Yule, (1). The Note of Fréchet adds comments 
on the measures r and n, listing certain limitations which follow from 
theoretical developments. In this note we add some comments which fol- 
low from the Regression Model, (1). Some additional references are given. 

1. It seems that Karl Pearson first used the Normal Probability Sur- 
face for two variables as a model, and defined the correlation r as the 
product moment of deviations of the variables from their respective means, 
with the standard deviations as units of measure. The two straight lines 


located on the plane of (X ,, X.,) by the method of least squares were used 
by Yule as the corresponding regression model for representing r, (1), 


page 188. The product formula for the frequency model is r = ais Yule 

gives the equivalent formula r = Vite from the regression equation z, = 


| 


<2? 


2 

2. When the points in the z,-arrays spread widely about the corres- 
ponding points of the regression line, the value of r is relatively small, 
and estimates of z , are not materially improved by using linear regression. 
When points in the arrays cluster closely about the corresponding points 
of the regression line, the value of r is relatively large, and estimates 
of x, may be improved by using linear regression. 


*r., where o ,, is the standard deviation of the linear estimates. 


3. The correlation ratio is not less than rfor a set of values(X ,, X ,) 
because the value of 7 is determined by the manner in which values within 
an array vary about the mean of the array, and the dispersion about the 
mean is less than that for any other point of an array. The values of r 
and 7 are equal for a given set (X ,,X,) if the means of arrays fall on the 
regression line. The values of r and y are abstract numbers calculated 
without regard to the manner in which values of X , and X, may be paired, 
Consequently they give no information regarding possible causes of cor- 
respondence. The worker must plan his problem so that some reasonable 
meaning may be assigned to the results. 
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4. Values of correlation may be relatively large in some cases where 
there is no apparent cause for correspondence. H. L. Rietz illustrated the 
extreme case where 7 is greater than r by using the function y = sing, 
O<a<2z.Here r=0 and 7 = 1. The comparison illustrates a case of non- 
linear regression. 

5. Yule, (1), page 195, refers to correlation as a measure of the 
strength of the influence of one variable upon another. The measure does 
not identify the influence. 

6. R.A. Fisher, (2), page 190, states the case as follows: “The cor- 
relation between A and B measures, on a conventional scale, the impor- 
tance of the factors which (on a balance of like and unlike action) act 
alike in both A and B, as against the remaining factors which affect A 
and B independently.” 

The case for common factors may be illustrated by the following ex- 
periment. Toss 5 coins on the table and count the number of heads (X |). 
Leave the 3 coins nearest the edge of the table and toss the other 2 coins 
down with them. Count the number of heads (X,). Repeat the experiment 
to obtain a sample of pairs (X ,, A..). The two series of counts are corre- 
lated and we say the 3 coins common to the two counts of heads consti- 
tute the common factors. In performing the experiment there could be other 
factors such as imperfect coins, or bias in the manner of tossing them. 

7. Rietz, (3) pp. 80-81, gives a probability mode! for correlation in 
the form y)dzdy = y)drdy. When A(z, y) is a function of both 
x and y, the two variables are said to be correlated, 
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MISCELLANEOUS NOTES 
Edited by 
Charles K. Robbins 


Articles intended for this department should be sent to Charles K. Robbins, De- 
partment of Mathematics, Purdue University, Lafayette, Ind. 


A FURTHER GENERALIZATION OF NEUSTADT’S LAW 


George Bergman 


In his article “A Mathematical Munchausen” in the Nov.-Dec. °56 
issue of this publication, Paul Schillo presented, amidst a large quantity 
of satire, a sequence of numbers called “Neustadt’s Base” and a certain 
conjecture, “Neustadt’s Law”, concerning this sequence. The sequence 
was: .2, -—l, 2, —2, 1, —.2; and the “law” that if this series ,be multi- 
plied, term by term by any succession of arithmetic progressions, the re- 
sulting sequence will have zero sum. E. g.: 

2 2 -2 -.2 
xl x3 
6 
x5 xl 
1 6 -15 
x1.80x1.44 «1.08  x.36 


1.80-8.64 6.48 5.76 -5.40 


1.80 -8.64+ 6.48 +5.76-5.40+0 = 0 


In the Jan.-Feb. °58 issue, W. W. Funkenbusch stated the correction 
that the number of successive multiplications must be limited to 4, and 
made some observations concerning multiplication by progressions of 
higher degree. 

However, all this is merely one case of a much more general law, 
with numerous applications, namely: 

If P is a polynomial of degree less than n, 


n 
= 0 
i=0 
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The proof of this consists of three parts : 
a) A ae P of degree less than n can be written in the form P(z) - 


a( ,)+ Q(z) where Q is a polynomial of than n-1 
To ion this, let A equal the coefficient of 2”~' in P(z) (which may 


be equal to zero). _— ia =) is a polynomial inz of degree n-1, in which 


the coefficient of . Thus, if we subtract (n-1 AC os from 


P(z), we will get a polynomial of degree less than n-1. L.e., letting 
(n—1)!A = a, we have 


P(x) - ,) = 


or 
P(x) = al, *,)+@(z) (Q of degree n-1) 
QED. 
b) 

The only values of ¢ for which (,," ,) 4 0 within the indicated range 
are n and n-1l. 

Thus 


nr 
i=-0 


c) If for all polynomials Q(z) of degree less than n-1, 


n—1 
> EW = 0 
7=0 
then for all polynomials P(x) of degree less than n, 


nr 
”)P(i) = 0 
i=0 


To show this we begin by applying a, and then b to the expression given: 
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n n n 
= >) (-1)( (-1)(?)Q) = (~1)"( 
i=0 i=0 i=0 
Now, as Q is a polynomial of degree less than n-1, we can prove this 


last expression equal to zero if we can break it into expressions of the 
form given in our hypothesis. This can be done as follows: 


n n 
> piped = 
t=0 i=0 


n n 
t= i=0 


The final term of the first of these sums is zero, and first term of the 
second sum is zero. We can adjust the ranges accordingly. To make the 
second conform with the form in the hypothesis, we shall then substitute 
j = i-1. Q(j+)) can be expressed as a simple polynomial in j, of the same 
degree as Q, which we shall call Q’: 


n 
pew 
t=0 


n-1 n 


i=] 
j=0 


n—1 
= 0- = 0-0 = 0 
j=0 
This proves c), which, together with the self-evident case when n = 0, 
constitutes an inductive proof of the theorem. 
Now “Neustadt’s Base” merely consists of the numbers (-1)%( 5), 
disguised by division by 5. Since the product of less than five arithmetic 


progressions is a progression whose general term is a polynomial of de- 
gree less than five, the sum of the products of corresponding terms of 
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less than five arithmetic progressions, multiplied by the corresponding 
terms of Neustadt’s Base, will equal zero divided by five, or zero. (“I. 
Neustadt” did not believe in mathematical induction — butone can’t please 
everybody! The particular case of “Neustadt’s Base”, or any particular 
case, can be proved deductively by multiplying the numbers of this “base” 
successively by a, a+d, a+2d, ---; then by a’, a*+d’, +--+; then a%ja°+d’j-«-.) 

The applications, as I mentioned, are numerous. For instance, take 
the method commonly used for finding the degree of a progression {S_}. 


We find first the differences between successive terms, AS, and then the 
differences of the differences, A?S, and so forth, until we reach an n such 


that all the A“S equal zero: whereupon we can state that n-1 is the de- 
gree of the progression. This can be explained in terms of the above 


theorem when one considers the fact that the kth A"S is 


n 
> 
i=0 


Another use for it is in the summing of progressions whose general term 
is P(2)r?.* Since the numbers (-1)(?) are the coefficients of the r° in 
the expansion of (1-r)’, if one multiplies the series by (1-r)°*+! where d 
is the degree of P(), then the coefficient of r* becomes 


d+1 

Pea, 

i=0 
or 0, except for the first and last d terms, where the sum will be incom- 
plete for want of non-included terms. Thus, all but 2d terms will drop out, 


and dividing through by (1-r)*, one will have a closed expression for the 
sum. 
Ezample : Evaluate 


S = 14+ 49+ + (n—2) + (n—1) 299-2 4 
Solution: Multiply through by 1-3r+3r?-+r* to obtain 


(l-r) 8S = 1-740 + 
(1-r)8 
*as, for instance, Proposal 333. r= > 


Stuyvesant High School, New York, N. Y. 


A NEW LOOK AT AN OLD PROBLEM 


Charles E. Wingo Jr. 


I have always been deeply interested in mathematics, and for the 
last ten years mathematics has been a source of recreation and enter- 
tainment for me. It was thru this that the old problem 


y?+2=7 


was revived. It was over fifty years ago that I first heard of it. What has 
stimulated interest in this problem is the fact that it has proved to be as 
much a puzzleas a problem. 

Equations of the form 


z7+y=a, 
were pronounced insoluble by Prof. Quimby of Dartmouth, except by Des- 
cartes formula for solving bi-quadratics. 
Those who are further interested in some of the solutions to the 
problem 
e*+y=1l, y7+2=7 


are refered to the Am. Mathematical Monthly, Vol. V, pages 291-2-3-4-5 
and Vol. VI, pages 13-37, where will be found about eight solutions, 
some of which are quite complicated; bi-quadratics, factoring, and other 
methods. 


In the problem as treated here. an effortis made to develope a method 
of solution whereby all problems of the form 


can be easily solved, where z > y and 2, y, a and 6 are positive inte- 
gers, provided of course that integral solutions exist. This means that 
arbitrary values of a and 5 can not be taken, as 
y2+2=4. 
The method derived here is general in scope, subject to the limi- 
tation mentioned above. Consider the two equations 
z*+y=a 
and 
y*+2=b. 
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We have 

= a-b 
or 
(1) (x-y)2+y-1) = a-b 


In my study of this problem over the past few years, my interest centered 
on equation (1) above. For the solution, if any was to be found, would 
hinge on this equation. For it is obvious that (a—5) is the product of two 
factors (x-—y) and (x+y-1), and the key to the solution was to find which 
of these equate to (r-y) and which to (r+ y-1). 

At this stage the problem rested for some time. Then some time later 


it flashed to me that Vary +y and Vy tte +2 could be written 


Vy*+2 ~ vb =yt+d 


Va-yb = r-y+e-d 


Since z > y, we know that in taking the square roots of r?+y and y?+2, 
the amount that c adds to z is less than the amount that d adds to y. 
Therefore d is larger than c. So Va - Vb will be less than the true value 
of z-y by the amount of d-c. So the value of the factor of a—4 to equate 


to z-y will be the integral value next greater than Va-\/b and (r+ y-1) 
would be equal to the other factor. 

To show how the method is applied to problems we consider the 
following. 


z*+y=ll=a 
= 7=6 
a-b=4 


Factors of 4= 1-2-2 
Va = = 3.31 
Vb=\7 = 2.67 
Va-yb = 
The first integral factor larger than .65 is 1. So 
z-y=1 
(continued on page 288) 


PROBLEMS AND QUESTIONS 
Edited by 
Robert E, Horton 


Readers of this department are invited to submit for solution problems believed to be 
new and subject matter questions that may arise in study, in research, or in extra-aca- 
demic situations. Proposals should be accompanied by solutions, when available, and 
by any information that will assist the editor. Ordinarily, problems in well-known text- 
books should not be submitted. 

Solutions should be submitted on separate, signed sheets. Figures should be drawn 
in India ink and twice the size desired for reproduction. 

Send all communications for this department to Robert E. Horton, Los Angeles City 
College, 855 North Vermont Ave., Los Angeles 29, California 


PROPOSALS 


376. Proposed by C.W.Trigg, Los Angeles City College. 

Identify the following unorthodox definitions : (1) the earthy or stony 
substance in which an ore or other mineral is bedded, (2) fully sufficient, 
(3) the four brightest components of @ Orionis, (4) type of speech fanci- 
fully exaggerated, (5) oddity, (6) a Franciscan friar, (7) a deposit formed 
in a liquid vegetable extract, (8) a large molding of convex profile, (9) 
fancied, (10) the series of air bubbles made by the breath from an otter 
under water, (11) a drain to carry off filthy water, (12) one sixteenth of a 
fluid drachm, (13) scheme, (14) depicted as walking, (15) absent-minded, 
(16) seekers, (17) complete, (18) invalid, (19) an interpreter of music. 


377. Proposed by J.M.Gandhi, Thapar Polytechnic, Patiala, India. 
Prove that 
9 3n+1 3 3n+2 
i=0 
when n is of the form 2m, 4m+3 or 4m—-1 respectively. 
378. Proposed by Barney Bissinger, Lebanon Valley College, Pennsyl- 
vania. 
Let a,, a, be arbitrary and n(n+l)a,,, = n(n—l)a,-(n-2)a,_, for 


n>0. Find S=Sa,. 
0 n 
379. Proposed by J.M.Howell, Los Angeles City College. 


Two men fight a duel. They both fire at a given signal. If both are 
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alive, they fire again at a given signal, and repeat the process until at 
least one of them is dead. If the probability that A kills B on any round 
is p, and the probability that B kills A on any round is r, find the proba- 
bility that A is alive, B is alive or neither is alive after: 1) n rounds, 
2) after an infinite number of rounds. 


380. Proposed by Huseyin Demir, Kandilli, Eregli, Kdz., Turkey. 
Solve the system of equations 


2(2—a) + u(z+u) = 0 

yla—b) + uly+u) = 0 

a(y—c) + ulz+u) = 0 
where abc 40 and 


381. Proposed by George M. Bergman and Melvin Hochster, Stuyvesant 
High School, New York. 

Show that if the circumsphere and the insphere of a tetrahedron are 
concentric, four spheres can be drawn each tangent to each edge of the 
tetrahedron, extended if necessary. 


382. Proposed by C.N. Mills, Sioux Falls College, South Dakota. 
For the hypocycloid 2?/*+ y?/3 = a?/%, («, 8) is the center of curva- 
ture. Show that «+8 = (2 


SOLUTIONS 


Late Solutions 


350. J.L. Brown, Jr., Ordnance Research Laboratory, Pennsylvania State 
University; M.Morduchow, Polytechnic Institute of Brooklyn. 


The Wine Merchant 


355. [November 1958] Proposed by P.B.Jordain, New York, New York. 

A wine merchant had a small cask containing 100 gallons of wine. 
In order to make more money, he decided to replace each gallon he took 
out of the cask by a gallon of water. This he did n times. Finding that 
he was losing customers, he naively tried to undo what he had done by 
selling his watered wine from the same cask, but replacing each gallon 
sold of the mixture by a gallon of pure wine, feeling that by this method, 
at the end of an additional n gallons of watered wine sold, he would have 
a cask of full wine again. He was unfortunate enough in having hit upon 
a total number of gallons sold such that the wine contained in the cask 
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was at a minimum at the end of this operation. The question is, how many 
gallons of watered wine had he sold in all? Assume that the merchant 
dealt only in full gallon sales. 


Solution by David F. Finlayson, Worcester Polytechnic Institute, 
Massachusetts. After n sales the fraction of wine remaining would be 


given by (578 )". The gallons of wine remaining would be given by 


100( =) = aoa: The gallons of water left would be given by 100- 
(99)" 


Gon ™=" After repeating this process n more times by replacing the mix- 

ture sold by pure wine instead of water as in the first n sales, the gallons 
99 (99)" — 

of water left would be given by The wine in the 

cask would then be given by 


_( 99 
100 )"1100- =. 


or 


G = 100-—(99)"__, _(99)” 
= * G00) 


G = 100 -(99)"(100) + (99)? 100) 


Differentiating with respect to n and setting & = 0 gives 
n 


=~ log, 99+ log,,100+ + (100)~"log ,100(09)"(-2) = 0 
n 


log .99 -log,100 = 2( )"(log99 log, 100) 


1 = 2( 
100 


0.5 = (0.99)” 
n = 69 * 
2n = 138 Qn-1 = 137 


Since the first gallon sold was pure wine, the wine merchant sold 
2n-1 = 137 gallons of watered wine. 

Substituting the value of n in the expression for the number of gal- 
lons of wine left, we find that the wine merchant ends with 83.57 gallons 
of wine in the cask. 
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*It can be immediately seen that (.99)®° can not equal exactly 0.050 
since the last digit of .99 to any power can not be 0. Solving the equa- 
tion by means of a slide rule or logarithms produces an answer very close 
to 69. It is better, however, to square both sides of ihe equation and 
solve for 2n. By doing this an irrational answer between two succesive 
integers will be found and of course the even integer is the answer de- 
sired. The value for 2n is then found to be 138 and n must equal 69 and 
not 68 or 70. 

Also solved by Stephen A. Andrea, Oberlin, Ohio; Philip Fung, Idaho 
State College; Joseph D. E. Konhauser, Haller Raymond and Brown, Inc., 
State College, Pennsylvania; Sam Kravitz, East Cleveland, Ohio; Richard 
Mittleman, Los Angeles City College; Lawrence A. Ringenberg, Eastern 
Illinois University; C.W. Trigg, Los Angeles City College; DaleWoods, 
Idaho State College; and the proposer. 

Several solutions were received with the value of n incorrectly ap- 
proximated by 68 instead of 69. Ringenberg pointed out that this problem 
is a special case of problem E962 appearing in the March - 1952 issue of 
the AMERICAN MATHEMATICAL MONTHLY. 


A Reciprocal Sum 


357. [November 1958] Proposed by Joseph Andrushkiw, Seton Hall Uni- 
versity, New Jersey. 
Show that if d>0, 


> 1/(k?-2dk+2d) = k+3a/4d, O<k<1/2d? 
k=0 


Solution by Arne Pleijel, Trolhattan, Sweden. Let the proposal be 
restated in the form: 


k?-2dk+2d? 4d ~ 2d? 
k-0 


Then we have 


dz 
2_2dx+2d? 
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k+1 


1 


and thus, as Gada? is decreasing in z, we have 


i 37 
—< 


1 
sat), (k 


and so the proposition is proved. 
Also solved by the proposer. 


The Bouncing Ball 


358. [November 1958] Proposed by C.W. Trigg, Los Angeles City Col- 
lege. 

\ ball having fallen from rest a vertical distance 4, strikes a stone 
protruding from a wall and bounces off horizontally without spinning. If 
the distance of the stone from the ground is s and the coefficient of res- 
titution is e, show that 

a) the ball will strike the ground at a distance 2\/she from the foot of 
the wall; 

b) the inclination of the stone to the horizontal is arctan, ?. 

Solution by Philip Fung, Idaho State College. Let u, v be the veloci- 
ties of ball before and after impact and their components paralle! and per- 
pendicular to the surface of stone be u,, vu, and v,, vy respectively. Let 
the inclination of stone to the horizontal be 0, the time taken by the ball 
to strike the ground after impact with stone ¢, and the horizontal distance 
covered by the ball during this time be PF. 

Now, the velocity of ball just before impact with the stone is u = 


V2gh. So 
u, = uSin 6; 


uy = uCos 6. 


Since uw, is parallel to the stone surface and thus unaffected in magnitude 
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and direction by impact, therefore we have, 
=u, = uSin 6; 


y= eu, = eu Cos 0. 


But, the direction of v is horizontal. Thus, 
Tan @ = = e/tané. 
From this we have, 
6 = arctanye. 
Now, 
s = 1/2 gt? 


or, 


t = (2879 


R=vt=t 


And, 


Hence, by substituting the values of ¢, v_, Vy and then @ into the equa- 
tion above, we can thus solve that R = 2\/she. 


Also solved by William E. F. Apphun, St. John’s University, New 
York; Joseph M. C. Hamilton, Los Angeles City College; Joseph D. E. 
Konhauser, Haller, Raymond and Brown, Inc.; State College, Pennsyl- 
vania; Dale Woods, Idaho State College; and the proposer. 


= 
Sing 
h 
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The String-of-Pearls Polynomial 


359. [November 1958] Proposed by Norman Anning, Alhambra, California. 

Prove that the string-of-pearls polynomial z* +1 can be expressed in 
at least one non-trivial way as the sum of two squares, if k is any even 
positive integer different from 2, 4, 8, 16---. 

Solution by William M. Sanders, Mississippi Southern College. Let 
k be the n-th term of an arithmetic progression with first term 0; i.e., k = 
(n—1)d. Impose the condition that d be the maximum divisor of & such that 
n = 2p > 4 to obtain the progression 0, d, 2d, ---, (n—2)d, k. This is pos- 
sible since & is even and different from 2, 4, 8, 16 --». Now for k = 2m, 
d = 2q, and n = 2p, 


2p—-2 
(-1)°4 min(v, 2p—v-1) 
v=1 
2p—2 
(-1)"4 min(v, 2p—v—1) 4 1 


v=l 


p-l p-l 3 

— 

v=l 


v=] 


Also sulved by L. Carlitz, Duke University; Arne Pleije!, Trolhattan, 
Sweden; and the proposer. 


A Difference Equation 
360. [November 1958] Proposed by Chih-yi Wang, University of Minnesota. 


Considering the higher differences of (*) with respect to z show that 


Solution by Arthur E. Danese, Union College, Schanectady, New 
York. In Rodrigues’ formula for Jacobi polynomials [Szego, ORTHOG- 
ONAL POLYNOMIALS, New York, 1939, p. 66, 4.3.1]: 


2°n! dz 


replace z by ¢-1, n by n-2z, B by z, « by r-n+m, and evaluate at ¢ = | 


to obtain 
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which equals 


y) = (-1)"P 


The use of the explicit evaluation of the Jacobi polynomial : 


n 
v=0 


{ibid, p. 67, 4.3.2] yields 


v=0 


which is equivalent to A"("), 


Also solved by the proposer. 
A Variable Nine-Point Circle 


361. [November 1958] Proposed by N. A. Court, University of Oklahoma. 

A variable triangle inscribed in a rectangular hyperbola has a fixed 
vertex and the opposite side moves parallel to itself. Show that its vari- 
able nine-point circle passes through two fixed points. 

Solution by J. W. Clawson, Collegeville, Pennsylvania. It is well 
known that the orthocenter of a triangle inscribed in a rectangular hyper- 
bola lies on the hyperbola. Also that the center of the hyperbola lies on 
the nine-point circle, (C. Smith, Geometrical Conics, pages 165, 166.) 

Let A be the fixed vertex, BC the variable side which moves parallel 
to itself. Let the perpendicular from A to BC intersect the hyperbola at 
H. Then the fixed pointH is the ortho-center of the variable triangle ABC. 
Since the nine-point circle of a triangle passes through the mid-point of 
AH, the variable nine-point circle passes through this point and also 
through the center of the hyperbola. 

Also solved by Huseyin Demir, Kandilli, Eregli, Kad., Turkey; 
Sister M. Stephanie, Georgian Court College, New Jersey; and the pro- 
poser. 


284 
m, 
since 


1959) PROBLEMS AND QUESTIONS 


Comment on Problem 337. 


337. (March 1958 and November 1958] Proposed by Victor Thebault, 
Tennie, Sarthe, France. 
Comment by William E. F. Appuhn, St. John’s University, New York. 
In order to eliminate having sides less than or equal to zero, condi- 
tions should have been placed as the solution was developed: 
1, In equation (1) z = 2pq, the condition pg > 0. 
2. In equation (2) y = »?-g?, the condition p? > p?. 
3. In the equations for p and g we do have to place restrictions on 
the integers g and A. They are: 
(a) g#0, (b) A £0, (c) g? # h?, (d) gh > 0 if g? > A?, and (e) gh < 0 
if g? < h?. 

. In the final set of equations, we have to place the restriction 
gh> 0. 
The above conditions will be satisfied by gh > 0 and g? > A?. 

>. Noting that, by equations (1) and (3), z and z are both even hence 
y is also even, a much quicker and easier complete solution may 
be obtained as follows: 

Let (1) = 4uv, 6 = and hypotenuse e = = 2d? 
where u and » are different positive integers. In order that u2+? = d? let 
(2) u- v = |x? -y?|, 0. d = where z ¢ y are integers. Substi- 
tuting (2) in (1) to obtain the immediate solution: 

(3) a = Bry |x?-y?|, b = and c = y?)?. 


QUICKIES 


From time to time this department will publish problems which may be solved by 
laborious methods, but which with the proper insight may be disposed of with dispatch. 
Readers are urged to submit their favorite problems of this type, together with the elegant 
solution and the source, if known, 


Q 247. For what values of z is m?+n?-a?~b6? > (mn-ab)z where 0 <a<m 
and 0<b<n. [Submitted by M.S. Klamkin} 


Q248. Show that if & is any real number, the lines 24+k2r°y-6r7y?- 
key*®+y* = 0 cut 2?+y? = 1 into eight equal parts. [Submitted by Norman 
Anning] 


Q249. Show that for the smallest values of the angles, 2 arccos 1/3 = 
arc cos (1/3) = Il. [Submitted by C.W. Trigg] 


Q 250. Prove that 19'° can not be represented as the sum of a fourth 
power and a cube. [Submitted by D. L. Silverman] 
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Q251. In a rectangular coordinate system (z, y), y = f(x). If the coordi- 

nates are transformed into new ones, (z‘’,y’), by a rotation through a posi- 
2 

tive angle «, s and int f 24 and [Submitted b 

ive angle «, express an in terms of and [Submitted by 

Henry E. Fettis] 

Q 252. Find the sum of 


[Submitted by M. S. Klamkin] 


Answers 


‘(usm ) = § OM 
= 


zp x 


ip pus jo ul pesseidxe aq yoiym wo 


Z/e le ele (ap) = 
= 


hep 


*SOXB JO OY} JOPUN SI By} 
(@p/Ap)+] UB] UBD + | 


“(ET Pow) 2 = = 

= 619 = 6l MG “(EL POW) 2 SuryzAue oq Avw 
SNYL “(ET Pow) 6 40 = puv (g] pow) ZI 40 “g ‘1 = 
= puB = D ‘asvo oy} 

= ¢/(DZ-9) +g /D\ S09 10 Z/4 = |] S09 
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> 10,09 < JSNU OM 


299 9S “9 = = Gy “Y= OV = gp ey) 
9 S09 0 <9 SOD UWZ BABY OM "9 = 
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(A New Look at an Old Problem — Continued from page 276.) 
z+y-l=4 
22=6 
y=2 
Next take: 
z*+y=1052=a 
y*+2= 816=5 
a-b = 236 
Factors of 236 are 2-2-59 = 4x59. 
Va = /1052 = 32.43 
Vb = = 28.55 
Va-\b = 3.88 
The integral factor which is next and larger than 3.88 is 4. Then 
z-y=4 
z+y-1=59 
22=64 = 32 
2y=56 y= 28 


For a further consideration of the problem, take the case where 
a= b, Then: 
z*+y=a 


y*+2=a 


=a 
=0 
Since z and y are both positive integers, (x+y-1) is also a positive inte- 
ger. Therefore (r-y) is equal to zero. Therefore z = y. So the equation 
above becomes : 
a(z+1) = 


It is seen from the equation above that a is the product of two consecu- 
tive integers. Hence a has such values as 1x2, 2x3, 3x4 ete. 


Amelia, Va. 
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